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The qq semirelativistic interaction in the Wilson loop approach 
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The complete qq semirelativistic interaction is obtained as a gauge-invariant function of the Wilson loop and 
its functional derivatives. The approach is suitable for analytic evaluations as well as for lattice calculations. 
Here we consider three different models for the Wilson loop non-perturbative behaviour and discuss the related 
semirelativistic dynamics. 



1. INTRODUCTION 

The derivation of the qq potential is a long- 
standing problem; in Ref. |^,|^ the static and 
the spin-dependent potentials were first obtained. 
Recently, working with a path integral represen- 
tation for the quark propagator in the external 
field, the complete qq potential was ob- 

tained as a gauge- invariant function of the Wilson 
loop and its functional derivatives H: 



dtVq^^i\og{w{r)) 
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where 



ig f dx^'A^{x) 



W{T) = P exp 



( ) is the normalized average over the gauge fields 
Ajj_ and 

{{f{A))) ^ {fiA)WiT))/{fiA)Wir)). 

The closed loop F is defined by the quark (anti- 
quark) trajectories Zi{t) (z2(i)) running from yi 
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to Xi (x2 to y2) along with two straight-lines 
at fixed time connecting yi to y2 and xi to X2. 
Moreover we have 



9 



(3) 



(((F^.(zi)F,,(z2)» - {(F,A^,))){{F^,iz,)))) 
= -^5^^3^((f..(.i))>, (4) 

where SS'^''(zi) = dz^Sz':^ — dz'^dz'^ . As a conse- 
quence we have that the actual form of the poten- 
tial in (|l|) is completely known once the Wilson 
loop behaviour is given. 

The semirelativistic qq potential can be written 

as 



(5) 



Vqq = Vb + VVD + VSD, 

with 



i\og{W{T)) = dt Vo{r{t)) + WD(r(i)) , 
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where r{t) = zi{t) — Z2{t) and the symbol { }w 
stands for the Weyl ordering prescription. The 
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spin dependent interaction appearing in (|l|) has a 
transparent physical meaning the first term being 
the magnetic interaction, the second the Thomas 
precession, the third the Darwin term and the last 
one the spin-spin interaction. However, it can be 
rewritten in the usual Eichten-Feinberg form |^] : 
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Si • S2Vi{r). 
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It is clear that all the dynamics is contained in 
the V^{r) of Eqs.(|)-(0) and these are unambigu- 
ous functions of the Wilson loop (cf. (P). Sim- 
ilar expressions for the Vi in terms of insertions 
on the static Wilson loop were obtained in Q]. 
These expressions are suitable for lattice evalu- 
ation and were used in [|| to obtain lattice pre- 
dictions for the semirelativistic interaction. How- 
ever, this is also the ideal framework to formulate 
hypothesis on the Wilson loop behaviour (and so 
on the confinement mechanism) to be checked on 
the lattice and on the experimental data. In the 
following we use three different models for the 
QCD vacuum to obtain an analytic behaviour of 
the Wilson loop in the non-perturbative region 
and to predict the semirelativistic quark dynam- 
ics. Comparison among the results allows us to 
get some insight in the mechanism of confinement 



2. THE CONFINING SEMIRELATIVIS- 
TIC DYNAMICS 

2.1. Minimal area lav^r model (MAL) 

In Ref. {W{T)) was approximated by the 
sum of a perturbative part given at the leading 
order by the gluon propagator D^^ and a non- 
perturbative part given by the value of the min- 
imal area of the deformed Wilson loop of fixed 



contour F plus a perimeter contribution V: 
i\og{W{T)) = -^g'^ dx^ j dx^2 



n 
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Denoting by u'^ = u^{s,t) the equation of 

any surface with contour F (s € [0,1], t G 

[ii,if], uO(s,i) = t, u(l,t) = zi(t), u(0,t) = 
Z2 (t) ) we have: 



ftp pi 
mm / dt ds 



\ dt ds 



dufj 
~dt~dt 



ds ds 



which coincides with the Nambu-Goto action. 
Up to the order 1 /m? the minimal surface can be 
identified exactly with the surface spanned by the 
straight-line joining (t,zi(i)) to {t,Z2{t)). From 
Eq. (^) Vqq is obtained [^J|] . In particular for the 
static potential we have 



Vo^~~+ar + C. 
S r 

See Tab. 1 for the complete results. 



(9) 



2.2. Stochastic vacuum model 

Using the non-Abelian Stokes theorem and the 
cumulant expansion it is possible to write 

{W{T)) 



P exp ( ig j dSf,^{u)F^„{u,xo) 



dS, 



x{Ff,iuiiui,xo) . . . Ff,.^^{uj,XQ))c 



(10) 



with ( )cum defined in terms of average values 
over the gauge fields ( ) and F^„{u, xq) being the 
path-ordered product of the field strength ten- 
sor F^y{u) times two Schwinger strings connect- 
ing the point u with an arbitrary reference point 
xo on the surface S appearing in the non-Abelian 
Stokes theorem. 

Equation ([l^) is exact. The first cumulant van- 
ishes trivially. The second cumulant gives the 
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first non-zero contribution to the cluster expan- 
sion. In the SVM |^] one assumes that in the 
context of heavy quark bound states higher cu- 
mulants can be neglected and the second cumu- 
lant dominates the cluster expansion, i. e., that 
the vacuum fluctuations are of a Gaussian type: 



log( W^(r)) = - Y ^ dS^. (u) ^ dSxp iv) 



with 



x{Fi_ii,{u,Xq)F\p{v,Xo))c 
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- — {{u - v)x5yp - (u - v)p5i,x) 
du,, 



_ g2 {Tr Fp,{0)Fp,{Q)) 



(12) 



36 D{0) + Di{0) 

Eqs. and (|l|) define the SVM for heavy 

quarks. The correlator functions D and Di are 
unknown. The perturbative part of Di, which is 
expected to be dominant in the short-range be- 
haviour, can be obtained by means of the stan- 
dard perturbation theory: 

= ^^'^^ J- 4- higher orders. (13) 

3tt x'^ 

One of the main features of SVM is to insert the 
information coming from the lattice inside the 
analytic model. Indeed the non-perturbative be- 
haviour of D and Di was evaluated on the lattice 
(see §) 

/3DLR(2.2) ^ de-*l"l, fiD\^(x'^) = die-*il"l (14) 

with 5 = {l± 0.1) GcV, d = 0.073 GeV^ 5i = 
(1 ± 0.1) GeV, di = 0.0254GeV^. From Eqs. 
(^l])-(|l^ all the potentials follow directly. In par- 
ticular the static potential is given by 



dr 



d\{r - A) 



xLI(t2 + 



j\\^D,{T^ + \^)Y (15) 



In the r — > limit (|T^) reproduces one gluon ex- 
change while in the limit r — > cxd 



/i-OO nOQ 
dr I dX D{t^ + A^) , 
-oo Ja 

/+OC poo 
dr dXX Dir^ + X^), 
-oo Ja 

/-foe pOC 
dr dXX Diir^+X"^). 
-oo Ja 

It is clear that (|l5| ) contains the MAL result 
for the static potential upon identification of the 
string tension and of the constants. For the other 
potentials see Ref. (complete form) and Tab. 
1 (r — > oo limit). 

2.3. Dual QCD 

The duality assumption that the long distance 
physics of a Yang-Mills theory depending upon 
strong coupled gauge potentials Ap is the same 
as the long distance physics of the dual the- 
ory describing the interactions of weakly coupled 



dual potentials Cp 
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CpXa/^ and monopole 



fields B^ = Y^ B^Xa/2, forms the basis of DQCD 

1^. The model is constructed as a concrete re- 
alization of the Mandelstam-t'Hooft dual super- 
conductor mechanism of confinement. Since the 
main interest is solving such a theory in the long- 
distance regime, the Lagrangian Ccs is explicitly 
constructed as the minimal dual gauge invariant 
extension of a quadratic Lagrangian with the fur- 
ther requisite to give a mass to the dual gluons 
(and to the monopole fields) via a spontaneous 
symmetry breaking of the dual gauge group. 

We denote by (Wcff(r)) the average over the 
fields of the Wilson loop of the dual theory M : 



(w^eff(r)) 



JVC VB ^ / ''^[•^""('^m^'+'^gf] 



JVCVB e 



i J dx[C,it{G%=0)+CGF] 
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£gf is a gauge fixing term and the effective dual 
Lagrangian in presence of quarks is given by 



U{Bi) is tlie Higgs potential with a minimum at 
a non-zero value Bqi = BqXt, B02 = and 
Bo3 = BqX2. Moreover we have taken Bi = B2 = 
B, the dual potential proportional to the hyper- 
charge matrix = C^Y and 



df_,Bi + ie[C^,Bi 



2tt 



gstiuap j ds j dr 



ds dr 
xd{x - y{s,T)), 



where j/(s, r) is a world sheet with boundary T 
swept out by the Dirac string. 

The assumption that the dual theory describes 
the long distance qq interaction in QCD then 
takes the form: 

{W{T)) = (Wcff(r)) , for large loops T. (16) 

Large loop means that the size R of the loop 
is large compared to the inverse mass {M^^ ~ 
(600 MeV)~^) of the Higgs particle (monopole 
field). Furthermore, since the dual theory is 
weakly coupled at large distances, we can eval- 
uate (Woff(r)) via a semiclassical expansion to 
which the classical configuration of the dual po- 
tentials and monopoles gives the leading contri- 
bution. In the leading classical approximation 

zlog(Weff(r)) = - j dx C,s{G%) , (17) 

with Cce{G^v) evaluated at the solution of the 
classical equations of motion. For the static po- 
tential we have 



This results directly from the classical solution 
to the nonlinear equations obtained from Ccs 
and reproduces the one gluon exchange in the 




limit r — > and the string term in the limit 
r — + c». The coefficient of the exponent can be 
actually given in terms of the dual gluon mass 
M = Qg'^Bl j-^- An interpolation of the nu- 
merical results for the potentials can be found in 
|^,||. In Tab.l we report the potentials in r cx) 
limit. 



3. FLUX TUBE STRUCTURE AND 
CONFINEMENT MECHANISM 

To discuss the semirelativistic non-perturbati- 
ve dynamics is convenient to study the large r 
limit in the three above models as reported in 
Tab.l. The MAL result is the realization of the 
intuitive Buchmiiller's picture of zero magnetic 
field in the flux tube comoving system: dl^/dr 
which is given by the magnetic interaction is 
zero in the non-perturbative region; similarly the 
velocity dependent potentials come from 

the consideration of the flux-tube energy. From 
Tab.l it is apparent that at the leading order 
in the long-range limit, neglecting exponentially 
falling off terms, the SVM contains exactly the 
MAL model results. In the spin dependent sector 
of the potential, both the SVM and DQCD not 
only reproduce the long-range behaviour given 
by the area law, but also give 1/r corrections 
to dVi/dr and dV2/dr. These corrections are 
equal in both models and very near to the ab- 
solute value of the constant term in the static 
potential (the SVM also supplies for the explica- 
tion of this fact) . This perfect agreement is abso- 
lutely not trivial and seems to be very meaning- 
ful, since it arises from two very different models 
in a region of distances in which the physics can- 
not be described by the area law alone. It is now 
clear that the vanishing of the magnetic part in 
the non-perturbative region takes place only at 
the leading order in the long-range limit. There- 
fore, working in a Bethe-Salpeter context, there 
is no need to assume an effective pure convolu- 
tion kernel which is a Lorentz scalar. Velocity 
dependent contributions to the quark-antiquark 
potential are important. In fact the string be- 
haviour of the non-perturbative interaction shows 
up when we consider the velocity dependent part 
of the potential and this is also what the data 
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Table 1 

Complete MAL potential and long distance SVM and DQCD potentials; from parameterization (^4|) we 
have: as = fl, C2 = fl; ^ = ^;^2 = ^-,£2 ^ see the values of d and 6 given in the text. The 
dual gluon mass — j^; typical values are a — 0.18; ag — 0.35. 
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-0.118crr + 0.275^5^ 
-0.177crr + 0.25375^ 



require. The velocity dependent structure which 
arises from the DQCD model differs slightly in 
the coefficients with respect to the area law be- 
haviour. In conclusion SVM and DQCD repro- 
duce the flux tube distribution measured on the 
lattice and the general features coming from the 
area law. Both give analytical expressions for 
the Wilson loop which describe the evolving be- 
haviour of {W(r)) from the short to the long dis- 
tances but not all predictions are equal in the 
two models in the intermediate distances region, 
in particular in the velocity dependent sector of 
the potential, and also in the spin-spin interac- 
tion as well as in AI4. The up to now available 
lattice data ||] confirm the MAL results but are 
not sufficiently accurate to discriminate between 
SVM and DQCD models. 
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